Spin Susceptibility of a 2D Electron System in GaAs towards the Weak Interaction 

Region 
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We determine the spin susceptibility x in the weak interaction regime of a tunable, high quality, 
two-dimensional electron system in a GaAs/AlGaAs heterostructure. The band structure effects, 
modifying mass and g-factor, are carefully taken into accounts since they become appreciable for 
the large electron densities of the weak interaction regime. When properly normalized, x decreases 
monotonically from 3 to 1.1 with increasing density over our experimental range from 0.1 to 4 x 
10^^ cm~^. In the high density limit, x tends correctly towards x ^ 1 ^-^d compare well with recent 
theory. 

PACS numbers; 71.18.+y, 73.40.-c, 73.43. Qt 



I. INTRODUCTION 

The recent interests in a possible correlated many par- 
ticle state Q, 01 very low electron density of a two 
dimensional electron system (2DES) have sparked a series 
of experimental investigations in different materials. 
0,0,111 In these experiments, the spin susceptibility, 
which is proportional to the product of effective mass, 
m*, and g-factor, g*, usually deviates from its nonin- 
teracting "bare" band edge value, tooSo, and increases 
drastically with decreasing electron density. The focus of 
these experiments is on the low density region, whereas 
X has not been pursued to high densities. While we ex- 
pect m*g* to approach mo.go in the high density limit, 
the explicit behavior in the intermediate regime remains 
experimentally unresolved. 

The density dependence of x is usually expressed in 
terms of the dimensionless parameter, Tg. It is defined as 
the average e-e distance measured in units of the Bohr 
radius, which is also the ratio of Coulomb energy to 
Fermi energy. In 2D, cx l/\/n, the intermediate den- 
sity regime resides at rg ~ 1 and the high density regime 
at Vs — > 0. High r^ measurements of m* g* from differ- 
ent materials generally differ from each other and from 
theory. Most early calculations |^ fiol Hit ll^ produced 
the sharp rise of x a-t low densities, but differed quanti- 
tatively from experiments. Of those, data on very thin 
AlAsQI quantum wells came closest to theory, probably 
since the latter omitted finite thickness effects. Indeed, 
recent calculation 0, IT^ ITHI | have shown that most of 
the discrepancies arise from the finite thickness of the 
2DES and the resulting modification of the Coulomb in- 
teraction. These calculations for x extend far into the 
weak interaction regime, yet there are presently no data 
available to test their validity. An additional incentive 
for such measurement is the fact that an empirical rela- 
tion for X, observed in the low-density regime 7], when 
extrapolated to high density would imply vanishing x, 
whereas one expects xl"' ^ oo) = 1. 



In this paper we provide high quality data on x through 
the intermediate interaction regime and towards small rs 
for comparison with recent, realistic theoretical calcula- 
tions. Once the data are corrected for the energy depen- 
dence of band mass and g-factor, theory and experiment 
reach a remarkable level of accord and the data show the 
correct limiting behavior for large densities. 



II. DEVICE STRUCTURE AND THE 
TILTED-FIELD MEASURING TECHNIQUES 

Our primary sample is a heterojunction-insulated gate 
field effect transistor (HIGFET). It consists of a 5/iTO 
buffer GaAs layer grown by molecular beam epitaxy 
(MBE) on a (001) GaAs substrate. A subsequently 
grown, 200 period superlattice of lOnm AZo.ssGao.er^s 
and 3nm GaAs acts as a diffusion barrier. It is followed 
by bf-im of undoped GaAs, which functions as the chan- 
nel. This layer is covered by 5nm of AlAs and 4/im 
of ^^o.33Gao.g7^s, being the dielectric material, topped 
by 25nm of heavily doped GaAs, which serves as the 
gate. The specimen is processed into a 600/zm square 
mesa. We use standard photolithography to define fif- 
teen Ni-Ge-Au contact pads. Our specimen has mobil- 
ities as high as 1 X l{f cm? /V s. The density, n, can 
be tuned by varing the gate voltage, Vg. For this sam- 
ple, n(10i°CTO-2) = 0.15 -h 16 X Vg[V). Such a HIGFET 
allows us to measure a wide range of densities in a sin- 
gle specimen, eliminating sample to sample parameter 
fiuctuations. For our sample, which has a density range 
spanning from 0.1 x 10^^ to 4 x lO^^cm^^, the rs value 
can be tuned from 0.9 to 6, covering a wide range of e-e 
interaction strength. 

All measurements were performed in a dilution refrig- 
erator with a base temperature of ~ iQmK . The sample 
was mounted on a rotating sample holder to control sep- 
arately the total magnetic field, Btot, and the field, B±, 
perpendicular to the 2DES. We used conventional low 
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FIG. 1: A typical gate sweep trace showing the coincidence 
condition for i = 1/2, at B± = 0.07T and tilt angle=86.5°. 
The arrow indicates the density for coincidence to be n = 
3.3 X lO^^cm"^. inset: schematic diagram for the i = 1/2 
coincidence, where spin splittings from the same and neigh- 
boring Landau levels are evenly spaced. 



frequency ac lock-in techniques with excitation currents 
ranging from InA to lOOnA, chosen to avoid sample heat- 
ing. 

We determine the spin susceptibility, x, by the tilted 
field methodH^, in the context of GaAs HIGFETs 7j. 
We briefly review the method, referring to Ref. Jj for de- 
tails : For a 2D system in a magnetic field, the spac- 
ing between Landau levels, Hujc = eHB±/m* , depends 
on the perpendicular field B±, while the Zeeman energy 
AEz — g*f^BBtot, depends on the total magnetic field, 
Bfot- When \g* ^J,BBtot\ = i^c, where the coincidence 
index, i, is an integer or a half-integer j^^i the Landau 
levels of the spin-up and spin-down electrons form a lad- 
der of evenly spaced energy levels. A schematic diagram 
of evenly spaced energy levels is shown for i = 1/2 as an 
insert to figure 1. 

The coincidence condition has a clear signature in the 
oscillation amplitudes of the magnetoresistance. In a 
density sweep at fixed 9 and B, and assuming m* and 
g* not to depend on density, such evenly spaced energy 
levels reveal themselves as uniform oscillation of uniform 
amplitudes. However, in general, m* and g* depend on 
electron density and therefore such equal amplitude os- 
cillations only appear at a particular density, n. Only 
at this particular n the value of m*g* = 2imoB±/ Btot- 
As an example, fig. 1 shows such a sweep for the i=l/2 
coincidence. 



III. RESULTS 



Fig. 2 shows the so determined n-dependent values of 
m*g*. Our data cover low values and, hence, persist 
into the low interaction regime, whereas earlier data by 
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FIG. 2: Spin susceptibility x versus electron density. Black 
squares are measurements for i = 1/2 in perpendicular field, 
Bj_ = O.IT; grey circles are data for i = l/2 and Bj_ = 0.07T. 
The error bars in x represents the uncertainties in angle. Er- 
rors in density are smaller than the symbol size. Open sym- 
bols represent data of same-shaped filled symbols after cor- 
rected for band structure effects. The extrapolation of x to 
Btot ~+ is plotted as a dashed line(see text). After this cor- 
rection the right scale becomes effectively x = m* g* /rribgb, 
which is used in fig. 4. Low density data from ref.Q are 
plotted for comparison. 



Zhu et al. (also shown) are limited to the high interaction 
regime(rs > 3). Our data were taken in two different 
perpendicular fields B^ = O.IT and B^ = 0.07T, both 
using the coincidence index i = 1/2. The error bars in 
the y-direction represent the uncertainties in angle mea- 
surements. For comparison, we also show the low-density 
data for the i = 1/2 taken at B±^ = 0.07T from Ref.jJ. 
The small difference between our B± — 0.07T data and 
theirs is most likely due to a slightly thinner wave func- 
tion thickness owing to the extra AlAs layer in our sam- 
ple(sce above). 

To put the main features of our data in perspective, 
we first follow the conventional way of normalizing the x 
data simply to the literature band edge values of r7io5o- 
We provide this scale on the right side of Fig. 2 Clearly the 
spin susceptibility decreases monotonically with increas- 
ing density from 3.4 to 0.88 for electron densities from 
0.1 X 10" to 4 X 10"cm"2 Qn this log-log plot, the data 
from each B± follow a power law behavior with a similar 
slope as observed earlier in the low density regime 0]. 
Obviously the data drop below x = 1 for high densities, 
which is unphysical and contradicts all theoretical predic- 
tions. However, bandstructure effects, which modify the 
"bare" mass and g-factor of the carriers, alter the power- 
law and prevent this from occurring. While such effects 
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FIG. 3: Data and fit for coefficient got from literature 
values [2^ I25I (see text and equation (1)). 

are negligible in the low density regime, they are appre- 
ciable at higher densities and need to be corrected for. 
In particular, in GaAs, the band mass mf, increases with 
increasing density, while the magnitude of the g-factor, 
\gb\, decreases, with the latter being dominant. 

IV. CORRECTIONS FOR BAND STRUCTURE 
EFFECTS 

The density dependence of the electron band mass in 
a 2DES in GaAs due to band nonparabolicity has been 
widely studied and can be readily calculated [21| . For our 
density range 0.1 x 10^^ to 4x lO^^cm"^, the enhancement 
of rrib grows from 0.3% to 5%. The density dependence of 
the g-factor is more complex. In a 2DES, this dependence 
was first observed by Stein et al. using Electron Spin 
Resonance (ESR) and was later studied systematically 
via ESR by Dobers et al.|23|. They deduced an empiri- 
cal formula gb{B, N) = gto-c{N + l/2)B, where N is the 
Landau level index, B is the magnetic field, and gbo and 
c are two coefficients that may change with electron den- 
sity. The observation was explained theoretically [^l ^ 
a combination of nonparabolicity of the band structure 
and an increase of the penetration of the electron wave- 
function into the AlGaAs layer. Since the g-factor of 
AlxGai-xAs changes considerably with Aluminum con- 
centration and has a different sign from go in GaAs, the 
magnitude of the g-factor shrinks with increasing pene- 
tration of the electron wavefunction. In order to correct 
our data, we use the existing measurements 
of gbo and c values from samples with electron densities 
from n = 1.5 X 10^^ to n = 3.3 x lO^^cm"^. Adding 
gbo = —0.44 for n — > 0, we interpolate and achieve the 
phenomenological dependence, 

gbo(n) = -0.44 - 2.75 X lO^^'n + 4.98 x IQ-^i^, (1) 

where n is in 10^°/cm~^(see fig. 3). The coefficient c 
is assumed constant, c=0.012, since it hardly varies in 
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FIG. 4: Normalized spin susceptibility enhancement x versus 
Vs. The open symbols are data corrected for band structure 
effects due to nonparabolicity. The extrapolation to Btot = 
is shown as a dashed line as in fig. 2. rrib and gb denote the bare 
m and g values after corrected for band structure effect. The 
solid curves of different grey scale are theoretical calculations 
from ref . 0, [ij, [ij . The topmost curve is for a zero-thickness 
2DES, whereas all other solid curves are for finite thickness 
GaAs systems. 

Ref. Il^l ■ While Dobers et al. used an Al concentration of 
~ 35%, the interface in our sample consists of 100% Al. 
From band structure calculations we can estimate that 
such a 100% Al concentration at the interface reduces the 
g-factor renormalization by ~ 1/2 as compared to ^ 35% 
Al. When combined with a smaller penetration into the 
AlAs layer, the final enhancement m* g* / rnhgt is affected 
by less than 4% at our highest Vg and has no observable 
effect for > 1.5. A similar estimate yields a varia- 
tion of c by less than 1% from Alo ^^GaQ Q^As/GaAs to 
AlAs/GaAs. These calculations are approximate. How- 
ever, we deem them sufficient since their impact on our 
data resides within the experimental error bar. 

The open symbols in fig. 2 represent all our x data cor- 
rected for the density dependence of mass and g-factor. 
As required, all data remain above x = 1 ^^id they seem 
to be levelling off as n increases. Since we have been 
working with a small spin polarization (i = 1/2) our 
data can be assumed to be close to the unpolarized(i.e. 
Btot — > 0) case, which is the condition used in most theo- 
retical calculations and therefore of importance for com- 
parisons with experiment. Based on the i?tot-dependence 
of our low-density data we can actually provide an ex- 
trapolation scheme to Btot — *■ 0, which, since it is a small 
correction, should hold over the extended high density 
range provided here. 

We observe that the B± = O.IT data and the B± — 
0.07T data are slightly offset from each other in the ver- 
tical direction. In particular, at lower densities both fol- 
low the previously observed powerlaw, which had allowed 
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for a simple extrapolation scheme|3| to Btot 0. This 
extrapolation, based on different coincidence indices, i, 
can directly be translated to different B± and identical 
i = 1/2 as is the case in the present measurements (see 
footnote 0]). The result is a constant multiplicative 
factor, which appears as a small downward shift on our 
log-log plot of Fig. 2 is shown as a dashed line. While 
this extrapolation scheme is derived in the range of the 
powerlaw dependence between % and n of our data, we 
expect it to apply approximately also in the range of sat- 
urating X towards higher densities. 

An enhancement of m* and a deduced increase of g* 
by a high in-plane magnetic field has been reported by 
Tutuc et al.[23 According to their observations, the mass 
enhancement for our largest in-plane field of 3.9T, at den- 
sity 3.6 X lO^^CTO"^, is approximately 8%. Extrapolating 
their g* enhancement to this density and field, we find 
it to be less than ~ 2%. Our x enhancement from the 
B ^ extrapolation is ~ 9%, which is close to the over- 
all ~ 10% increase deduced from Ref . [2^ . This supports 
that our method, which extrapolates the spin suscepti- 
bility to vanishing B field, xb-^o — ^* g* /^nogo, is not 
sensitive to in-plane field effects on m* and g*. 

Fig 4 shows our final results as x versus , rather than 
versus density, for better comparison with theory. The 
symbols and lines refer back to fig. 2. The data points as 
well as the extrapolated curve show the correct limiting 
behavior of x ^ 1 as ^ as a consequence of the 
correction of the data for band structure effects. Without 
such correction, the x data would have fallen below x = 
1, which would be clearly unphysical. 

V. COMPARISON WITH THEORIES 

Finally, we can compare our data with theory. In the 
past two decades, there have been several calculations of 
the spin susceptibility x, but most of them evaluate x 
in the limit of zero polarization in an ideal 2DES. A fig- 
ure summarizing those theoretical results can be found in 
Ref. _9J . The calculations from different techniques are in 
rather good agreement for < 3 but deviate consider- 
ably for large . In the same regime, measurements of x 
in different 2DESs have also shown considerable discrep- 



ancies. A recent report of dePalo et al, attributes many 
of these discrepancies to the finite thickness of the real 
2DES and variations of this thickness between 2DESs 
from different material systems. For comparison with 
our data, extrapolated to B ^ (dashed line), we have 
included in fig. 4 the zero thickness result and the finite 
thickness result for GaAs from reflsj. The calculations 
from Zhang and Das Sarma 14J using an improved RPA 
and Dharma-wadanall5| using a CHNC method both re- 
side further below our data. 

As can be seen from fig. 4, our experimental data and 
the theoretical calculations for a realistic 2DES in GaAs 
show rather good agreement. For high rg values, this 
had been pointed out before [l^. Here we show that after 
appropriate correction for band structure effects, one also 
achieves good agreement in the intermediate interaction 
regime {vs ^ 1) and obtains the correct trend for Vg — > 0. 



VI. CONCLUSION 

In conclusion, our data on the spin susceptibility of a 
2DES in GaAs over a very wide of electron densities, cov- 
ering values from strong electron interactions (ts 6) 
through the intermediate interaction regime (r^ 1), 
and into the low interaction regime show the correct 
trends for the limiting behavior {r^ — > 0) of this quantity. 
Only when theory takes into account the finite thickness 
of the real 2DES and experiment is corrected for band 
structure effects, a rather good agreement between data 
and calculation is reached over the whole range of of 
this study. Our work highlights the fact that subtleties of 
the real 2DES need to be taken into account when com- 
parison with theoretical calculations are being made. 
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